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Abstract 

>■ We present an analysis of comet activity based on the Spitzer Space Telescope component 

"^ of the Survey of the Ensemble Physical Properties of Cometary Nuclei. We show that the 

QQ survey is well suited to measuring the activity of Jupiter-family comets at 3-7 AU from the 

cn Sun. Dust was detected in 33 of 89 targets (37 ± 6%), and we conclude that 21 comets 

"^ (24 ± 5%) have morphologies that suggest ongoing or recent cometary activity. Our dust 

i^s detections are sensitivity limited, therefore our measured activity rate is necessarily a lower 

^-H limit. All comets with small perihelion distances (g < 1.8 AU) are inactive in our survey, and 

>• the active comets in our sample are strongly biased to post-perihelion epochs. We introduce 

k> the quantity efp, intended to be a thermal emission counterpart to the often reported Afp, 

^ and find that the comets with large perihelion distances likely have greater dust production 

rates than other comets in our survey at 3-7 AU from the Sun, indicating a bias in the 

discovered Jupiter-family comet population. By examining the orbital history of our survey 

sample, we suggest that comets perturbed to smaller perihelion distances in the past 150 yr 

are more likely to be active, but more study on this effect is needed. 
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1. Introduction 

The Survey of the Ensemble Physical Properties of Conietary Nuclei (SEPPCoN) provides 
a data set with which we may examine comet activity at intermediate heliocentric distances, 
Th (here, 3-7 AU). The primary goal of SEPPCoN is to measure, in a consistent manner, the 
sizes and surface properties of a statistically significant number (^ 100) of Jupiter-family 
comet nuclei. SEPPCoN is comprised of two observational campaigns: a mid-infrared (mid- 



IR) component utilizing the imaging capabilities of the Spitzer Space Telescope (Werner et al. 



2004), and a visible light component using ground-based instrumentation. Early results from 



the Spitzer survey on Comets 22P/Kopff and 107P /Wilson-Harrington were presented by 
Groussin et al. (2009) and Licandro et al. ( 2009[ ), respectively. Fernandez et al. (2011), here- 
after Paper I presented the first survey results from SEPPCoN. In Paper I we measured 
and analyzed the thermal emission from 89 comet nuclei, and presented results on the cu- 
mulative size distribution and spectral properties of the Spitzer survey targets. The goal of 
the present paper is to examine the dust and activity of those same targets. 

Active comets produce a coma. A coma is a gravitationally unbound atmosphere, driven 
by the outflow of volatiles from the nucleus surface or sub-surface. As comets approach the 
Sun, their surfaces warm, eventually causing gases to be released from the nucleus due to 
sublimation of volatile ices. The heliocentric distance at which this occurs depends on the 
physical properties of the nucleus in question: shape, composition, and internal structure. 
Volatile driven mass loss due to insolation is typical of comet activity, i.e., we may not even 
consider an object to have cometary activity unless we observe a coma and/or tail generated 
by solar heating of volatiles (as opposed to impact driven mass loss). In addition to the 
gases, the coma typically includes escaping dust and/or icy grains lifted off the nucleus by 
the gas outflow. As a comet recedes from the Sun, the nucleus cools and activity may be 
quenched. 



Water ice is typically the primary driver of comet atmospheres inside of 3 AU ( Meech and 



Svoreh 2004). Water is by far the most abundant ice near the nucleus surface, as inferred 
from remote spectrophotometric compositional studies of comae. The next most abundant 



gases are CO2 and CO with relative fractions ^ 20% ( Bockelee-Morvan et al.||2004 Ootsubo 



et al.||2012 ). However, CO2 can also drive activity, even though it is generally less abundant 
than water, as shown by the Deep Impact flybys of Comets 103P/Hartley 2 (A'Hearn et al. 
2011) and 9P/Tempel 1 (Feaga et al. 2007). Outside of 3 AU, the relative contributions of 



water ice and more volatile ices to activity is less well understood. Activity has been observed 
in more than 80 comets with large perihelion distances {q > 5 AU), despite the inferred low 
surface temperatures at such great distances. Ices more volatile than water must play larger 
roles in driving activity for these distant comets than for comets closer to the Sun. Water 
sublimates at T ^ 170 K (which occurs in the solar system near rh ^ 3 AU), whereas CO2 
ice sublimates at T Ri 80 K {vh ^ 12 AU) and CO at T ^ 20 K [rt ^ 50 AU). Recent results 
from the Akari satellite show that the coma mixing ratio of CO2 to H2O is systematically 
larger for comets outside of 2.5 AU, indicating the diminishing role of water sublimation as 



heliocentric distance increases (Ootsubo et al. 2012). In addition to sublimation of ices, the 



crystallization of amorphous water ice could release trapped gases and drive activity at 120- 
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160 K ( Schmitt et al.||1989 Meech and Svoren||2004 ) , and has been proposed as the dominant 
driver of activity in Centaurs (Jewitt 2009). Whatever the mechanism that drives activity 



in comets with g > 5 AU, not all comets are active at such large heliocentric distances. By 
studying the activity of comets at intermediate heliocentric distances we may learn which 
properties of comet nuclei initiate activity as they approach the Sun, and quench or sustain 
activity as they recede from the Sun. 

Mid-IR broadband images of comets, such as those taken as part of our survey, are dom- 
inated by thermal emission from the nuclei and surrounding dust. With few exceptions, the 
gas is undetected or only forms a minor component of the emission. Because gas expansion 
is a critical component to comet activity, but remains undetected in most mid-IR images, in 
this work we infer activity from the presence and morphology of dust structures larger than 
the point source nucleus, rather than from the direct detection of gases. 

With SEPPCoN Spitzer images, we measure the dust activity of comet nuclei at 3- 
7 AU, distances where water ice sublimation is typically low. First, we review the Spitzer 
observations which were presented in detail in Paper I (^. Next, we present the morphology 



and the photometric properties of the dust detected in the survey, and assess the nature of the 
dust (Q. Then, we discuss the frequency of cometary activity versus heliocentric distance 
and other parameters, their implications on the structure and heating of comet nuclei, and 
analyze the color temperature of comet dust at 3-7 AU (Q. Finally, we summarize our 
findings (^. 



2. Observations and Reduction 

The purpose of the Spitzer component of the survey is to obtain a robust estimate of the 
size distribution of known Jupiter-family comet nuclei. To this end, 100 JFCs were selected 
that met the following criteria: 1) the ephemeris was constrained well enough such that the 
comet could be expected to lie within a 5' x 5' field of view (the footprint of Spitzer^s largest 
arrays); 2) the comet must have been observable by Spitzer and beyond ~ 4 AU from the 
Sun during Spitzer Cycle 3 (July 2006 to July 2007); and, 3) the nucleus must have been 
brighter than V = 24.0 mag to make ground-based optical observations feasible. Criterion 3 
requires an estimate of the nucleus radius (R), which we obtained from the compilation of 
Lamy et al. (2004). If no estimate existed, we used the following assumption: R = 1.0 km 
for comets with perihelion distances q < 2.0 AU, R = 1.5 km for 2.0 < q < 2.5 AU, and 
R = 2.0 km for q > 2.5 AU. The assumption attempts to account for the fact that it is 
increasingly difficult to discover smaller comets at larger perihelion distances. 

Two Spitzer instruments were well suited for the survey: the 24-pm (Ag// = 23.7 ]am, 
2.55 arcsec pixel"^, 5.4' x 5.4' field of view) camera of the Multiband Imaging Photometer 
for Spitzer (MIPS; Rieke et ar||2004 ), and the 16-iim and 22-]im Infrared Spectrograph (IRS; 
Houck et al. 2004) peak-up arrays (Ag// = 15.8 and 22.3 pm, 1.85 arcsec pixel"^, 0.9' x 1.4' 



field of view). Our choice of instrument was based on each comet's ephemeris uncertainties. 
If the 3a uncertainty was under 30" we selected the IRS peak-up arrays; if the uncertainty 
was between 30 and 200" we selected the 24-]am MIPS camera. The two IRS peak-up arrays 
were preferred because they allowed us to measure a color for each nucleus, which provides 
an indication of the effective temperature of the surface (required for the size estimate). 



The Spitzer spacecraft tracked each comet given the computed ephemerides, and observed 
each target twice. For the IRS, each target was first observed with the 16-pm peak-up array, 
followed immediately with the 22-pm peak-up array. IRS peak up background observations 
were obtained in the array that was not centered on the comet, although the background ob- 
servations were not always useful due to variable detector artifacts (especially latent charges), 
and nearby background sources (stars). For the MIPS, a duplicate (shadow) observation was 
executed 1 to 30 hr after the primary observation. The shadow observation was close enough 
in time to the primary such that the comet remained within the MIPS field of view but 
displaced from the original position. The goal signal-to-noise ratio on the total flux of each 
nucleus observation was 30, based on the radius estimates discussed above. In Table [T] we list 



the 89 comets identified in Paper I with their observing circumstances and relevant orbital 



parameters. Details on the general image reduction, and the identification of specific targets 



are presented in [Paper I[ Out of the 100 targeted comets, two were not in the field of view 
of their observations, six were in the field of view but were not detected, and three have an 
unknown status (i.e., the ephemeris is sufficiently uncertain that we cannot conclude if the 
comet was in the field of view) . 

3. Results 

3.1. Dust morphology 

There are three types of comet dust morphologies relevant to this paper: comae, tails, 
and trails. The coma begins at the surface of the nucleus. In the rest frame of the nucleus, 
entrained dust grains move away from the surface. At a telescope we observe a roughly 
elliptically shaped coma with a surface brightness distribution that decreases with distance 
from a central source. The presence of this extended dust coma is the best evidence for 
recent and ongoing activity in mid-IR observations. 

As the coma expands, the morphology becomes increasingly dependent on the dynamical 
and physical properties of the dust. Gas expansion accelerates dust grains from the nucleus 
surface and places each grain into a new orbit around the Sun. The grain's orbit depends on 
its velocity and response to solar radiation pressure. The radiation force is size dependent. In 
general, smaller grains feel a greater acceleration from radiation pressure. The trend reverses 



in the sub-micron range where grains are too small to efficiently absorb solar radiation (Burns 



et al. 1979). In studies of comet grain dynamics, the radiation force is commonly expressed 
as the unitless parameter /3 = Fr/Fg oc Qpr/a, where Fr is the force from solar radiation, Fg 
is the force from solar gravity, Qpr is the grain-dependent radiation pressure efficiency, and 
a is the grain radius. 

The nucleus may also experience non-gravitational forces. The comet nucleus does not 
feel an appreciable radiation force, but instead insolation drives mass loss that produces 



significant secular changes in the comet's orbit (Whipple 1950, Marsden et al. 1973). Alto- 
gether, the nucleus and the dust grains are in different heliocentric orbits, and as the dust 
coma expands it transforms into a dust tail. The presence of a dust tail may indicate recent 
activity, but it is not as conclusive as the presence of a dust coma. Ambiguity is present be- 
cause the tail, as projected on the sky, may be composed of intermediate-sized slow-moving 
grains that were ejected weeks prior to the observation, whereas smaller grains will have 
much shorter lifetimes. 



The largest grains are removed from the nucleus with the lowest ejection velocities, and 
weakly interact with solar radiation (i.e., they have small /3 values). The heliocentric orbits 
of the large grains are so similar to the nucleus that only months or years after ejection will 
their presence be apparent as a long, linear dust feature along or near the projected orbit of 



the nucleus (Sykes and Walker 1992). Such a linear structure may be interpreted as a dust 
trail, and its presence indicates activity on month-to-year timescales. 

For each image, we searched for the presence of dust using two methods: PSF matching 



and visual inspection. In Paper I[ we fitted the central source of each comet with a scaled 



PSF. Those comets with residual emission in excess of a smooth background were identified 
as potentially having dust. Visual inspection of each image served to verify the initial results 
from the PSF fitting, and to identify dust emission outside of the PSF fit radius (typically 5 to 
8 pixels). The images were inspected with a range of color scales and smoothing techniques 
to verify potential dust. Since there are two images of each comet, finding dust in both 
images gives credence to our interpretation, but is not strictly required as contamination 
from artifacts and background sources could obscure dust in one image but not the other. 
In practice, dust is easily found via visual inspection for most targets. Only in images of 
comet 50P/Arend were the PSF fit residuals suggestive of dust that could not be verified 
with visual inspection due to confusion with background sources (star streaks). 

If dust is detected, we describe it as a coma, tail, trail, or some combination of the three. 
To aid our identifications, we generated a set of zero-ejection-velocity syndynes (curves of 
constant (3 but variable emission time), calculated for /3 = 1, 10~^, 10~^, 10"'^, and 10~^, 



using the software of Kelley (2006) and Kelley et al. (2008). The survey images with dust 



detections are presented in Figs. [T] through [6} along with model syndynes. Our identifications 
are listed in Table |2] and are based on the following guidelines: 

1. Comets with bright dust surrounding the nucleus are described as having comae (e.g., 
32P, 74P, 213P). 

2. Comets with dust following the /3 > 10"'^ syndynes are described as having tails (e.g., 
159P, 173P). 

3. Comets with long linear dust features following the /3 < 10^'^ syndynes are described 
as having trails (e.g., 74P, 219P). We also labeled the dust as a trail if it leads the 
comet along the projected orbit (e.g., 22P, 144P). 

4. Comets with thin linear dust features that overlap multiple syndynes and the projected 
orbit of the comet are ambiguous, and are described with the label tail/trail (e.g., 56P, 
78P). 

5. Comets with broad, but short, dust detections along any syndyne are labeled as tails 
(e.g., 119P, P/2005 JDios). 

Out of 89 comets, 56 comets (63%) do not have clear morphological evidence for dust. 
The remaining 33 dust detections are described in Table |2] In addition, 2 comets are marked 
as tentative detections of dust: Comet 6P/d'Arrest may have a very faint trail along the pro- 
jected orbit of this comet, and Comet 50P/Arend has a slight surface brightness excess along 
the /3 = 0.1 syndyne (but there is also a nearby star, which complicates the interpretation). 
Both comets are shown in Fig. [Tj These cases are "tentative" as opposed to "ambiguous" 
because dust has not been definitively detected. 
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Figure 1: 5piteer/IRS 22-]jm peak- up images and model syndynes for comets with resolved dust (the 16-]jm 
images have similar morphologies, but lower signal-to- noise ratios). There are two panels for each comet: 
the lower panels are filled contour plots of Gaussian smoothed {<J = 1.5 pixel) versions of the upper panels. 
Also plotted in the lower panels are /3 = 0.1 (dotted lines) and /3 = 0.001 (dash-dotted lines) syndynes, 
and the projected orbit of the comet (solid lines). The nuclei are located at the intersections of the three 
lines. The right images are the 22-]jm data obtained while the 16-]am array was centered on the comet. In 
some cases, these bonus images contain dust (e.g., 22P). The top and bottom panels are plotted with the 
same reverse gray-scale data limits. The image orientations are indicated: Celestial North (N) and East (E), 
projected sunward direction (0), and projected velocity (v). Comets 6P/d'Arrest, 22P/Kopff, 32P/Comas 
Sola, 37P/Forbes, 48P/Johnson, and 50P/Arend are shown. 
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Figure 2: Same as Fig. [11 but for Comets 56P/Slaughter-Burnham, 69P/Taylor, 74P/Smirnova-Chernykh, 
78P/Gehrels 2, lOlP/Chernykh, and 118P/Shoemaker-Levy 4. 
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Figure 3: Same as Fig.[T] but for Comets 119P/Parker-Hartley, 121P/Shoemaker-Holt 2, 129P/Shoemaker- 
Levy 3, 152P/Helin-Lawrence, 159P/LONEOS, and 171P/Spahr. 
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Figure 4: Same as Fig.[T] but for Comets 173P/Mueller 5, 213P/2005 R2 (Van Ness), 246P/2004 F3 (NEAT), 
P/2004 V5-A (LINEAR-HiU), and P/2004 VRg (LONEOS). 




Figure 5: Same as Fig. Ill but for Spitzer /MIPS 24-]jm images. Most images are background subtracted 
with the second (shadow) observation of the same comet. In these cases, the grayscale data Umits are chosen 
to enhance the contrast of only one of the two comet images. Comets 16P/Brooks 2, 62P/Tsuchinshan 1, 
144P/Kushida, 219P/2002 LZn (LINEAR), 260P/2005 K3 (McNaught), and C/2005 W2 (Christensen) are 
shown. 
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Figure 6: Same as Fig. [s) but for Comets P/2004 Al (LONEOS), P/2004 H2 (Larsen), P/2005 JDios 
(Catalina-NEAT), P/2005 T5 (Broughton), P/2005 W3 (Kowalski), and P/2005 Y2 (McNaught). 
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3.2. Dust photometry 

The Spitzer images allow us to estimate each comet's dust mass and, for the IRS observa- 
tions, color temperature. For comets with a coma and/or tail, we measure the total thermal 
emission in 3-6 pixel radius apertures centered on the comet. The smaller aperture sizes 
are necessary in some images to reduce background contamination. When possible, we also 
examine an aperture offset at least 5 pixels from the nucleus in order to measure emission 
from tails and trails. The size and shape of the tail/trail aperture varies comet-by-comet 
and depends on the morphology of the dust. The background is estimated from areas of 
low contamination from dust and background sources. As the thermal emission from comets 
typically contains significant emission from the nucleus, we only use the nucleus subtracted 



images from [Paper I[ The fluxes have been aperture corrected and color corrected assuming 
the spectral shape of an isothermal blackbody sphere in local thermodynamic equilibrium 
(LTE), according to the methods prescribed by the MIPS and IRS instrument handbooks 



(Colbert 2010 Spitzer Science Center 2009). The dust fluxes are listed in Table |2| 

With the two IRS peak-up arrays, we measure the 16-22-pm color temperature of the 
dust in our survey, centered on the nucleus and in the tail/trail (if applicable) for all fluxes 
with signal-to-noise ratios > 5. Dust color temperatures are listed in Table [3] as a ratio 
with the temperature of an isothermal blackbody sphere in LTE at the same heliocentric 
distance {Tbb = 278 r^ ). The error- weighted mean of all color-temperature measurements 
is {T /Tbb) = 1.074 ± 0.006, of the center apertures is 1.079 ± 0.007, and of the tail/trail 
apertures is 1.065 ± 0.010. 

The color temperature of the dust will only serve as an approximation to the true physical 
temperature of the grains because physical temperatures depend on the grain sizes, shapes, 
and compositions, and comet dust structures are collections of grains with a wide range of 
properties. For example, dust grains with sizes < 1 \\.m have higher equilibrium temperatures 
than larger grains because the former have smaller emission cross sections at 10-20 pm, the 
wavelength regime at which the peak of the thermal emission occurs for a blackbody sphere 
in the inner solar system. Thus, their equilibrium temperatures are increased over that 
of a blackbody. These increases are composition dependent as carbonaceous grains, which 
have high absorption and emission efficiencies, will have different equilibrium temperatures 
than an equally-sized silicate grain, which is less absorptive at most optical and infrared 
wavelengths. For more discussion on these effects, and how they may constrain thermal 



emission spectra, see Wooden (2002). 

Despite the ambiguity in correlating color temperatures to coma grain properties, comet- 
to-comet differences indicate variations of those (unknown) dust parameters among the mem- 
bers of our survey. For example, the color temperature of 74P/Smirnova-Chernykh's coma 
is greater than lOlP-A/Chernykh's coma (1.16 ±0.02 versus 0.96 ±0.03, respectively). This 
difference indicates that 74P's coma is composed of different grains than lOlP's coma, and 
we can speculate that they are smaller and/or more carbon rich in the former. 

3.3. Survey sensitivity 
3.3.1. All dust 

The SEPPCoN survey was designed to detect comet nuclei with a signal-to- noise ratio of 
30. The exposure times are not based on potential dust emission but solely on our estimated 
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nucleus sizes and temperatures. To verify that we can draw meaningful results on the dust 
activity of Jupiter-family comets as a whole, we must estimate the survey's sensitivity to 
dust. Rather than using a theoretical estimate (i.e., observation planning tools), we prefer 
to measure the sensitivities from our observations, which will easily take into account the 
effects of background objects (stars, asteroids, etc.), variances in instrument calibration (e.g., 
flat-fielding), latent charge on the detector, etc. To measure the sensitivity of an image, we 
first filter the image to remove structures larger than 5x5 pixels (i.e., dust and the smoothly 
varying background). To be more specific, we applied a morphological gray closing operator, 
followed by the opening operator, and subtracted the result from the original image. The 
opening and closing operators effectively convolve the image with a 5 x 5 box, but instead 
of replacing each pixel with the average of the surrounding pixels, the closing and opening 
operators replace each pixel with the maximum and minimum pixel values within the box, 
respectively. Applying the closing and opening operators is similar to applying a median 
filter, except more of the small scale structure is lost in median filtering. Because small scale 
structure affects our abihty to detect dust, we prefer the closing and opening filters over 



the median. See Lea and Kellar (1989) and Appleton et al. (1993) for further discussion 



and other applications of morphological operators in astronomy. In the part of the co-added 
image where the integration time was the highest, we measure the standard deviation of the 
image after iteratively removing 3a outliers (i.e., the comet itself, and the central cores of 
stars). 

The image sensitivities are measured in units of MJy sr^^. However, dust surface bright- 
ness depends on heliocentric distance: a surface brightness of 1 MJy sr"^ at 3 AU implies 
less dust than 1 MJy sr~^ at 7 AU due to the different equilibrium temperatures. A more 
relevant representation of the sensitivity is needed. We have chosen to transform the mea- 
sured surface brightnesses into optical depths. To compute the image sensitivity to dust in 
terms of optical depth, we use the equation: 



CiT,)B,{T,)VA 

where a,- is the optical depth (or effective fill factor) of a la per pixel detection of dust, o"/^ is 
the measured sensitivity (la per pixel) in units of MJy sr~^, T^ is the effective temperature of 
the dust, C{Td) is the instrument color correction for a blackbody spectrum at a temperature 
Tfii B,^{Td) is the Planck function evaluated at temperature T^ in units of MJy sr~^, and A 
is the area of a fictitious dust aperture in units of square pixels. Our default aperture size 
in Table [2] is a 6 pixel radius circle. In ^3^, we computed a mean color temperature of 



Tc/Tbb = 1-08 in the center aperture. A priori, we do not know the color temperature of 
any specific comet. Rather than using one value for 70 comets and our measured values for 
those 10 comets in Table 3, we assume T^ ^ Tc = I.OSTbb = 300 r^ K for all comets. 

Histograms of the observed dust optical depths and the survey image sensitivities are 
presented in Fig. [7] for the IRS 22-pm and MIPS 24-]im observations. Notice that the dust 
detections fall off at the same optical depths as our estimated 3a image sensitivities. This 
fall off strongly suggests our dust detections are sensitivity limited. 



13 



3.3.2. Trail dust 

In a survey of 34 comets taken with the MIPS instrument, 
that at least 80% of all short-period comets have dust trails. Since dust trails are long lived 
features (^3.1) we may expect a similar rate for the SEPPCoN targets, but instead we find 



Reach et al. (2007) found 



a rate of 10%. 

To better understand the dramatic difference in trail rates, we compared our survey 



target list to that in Reach et al. and found the following 12 comets in both surveys. 



• Comet 107P did not have a trail, or any other dust, in any observation. 

• Comet 62P had a trail in both surveys. 

• Comet 78P's dust was classified as a possible trail in SEPPCoN, and as having a trail 
by [Reach et alT 



• Comet 121P had a trail in SEPPCoN but only a tail in the [Reach et al.| survey. We 
suggest that this comet's trail only becomes apparent at large true anomalies (trail 
survey / = 31°, SEPPCoN / = 123°). 



Comets 32P, 69P, 123P, and 131P had trails in the Reach et al. survey, but they 



appear to be too faint to detect in the SEPPCoN images. In particular, 123P's SEP- 
PCoN images have several nearby point sources and background subtraction artifacts 



(Fernandez et al. 2011). 



Comets 48P and 129P had clear leading and following trails in the Reach et al. survey. 



The fields-of-view of the SEPPCoN observations are limited in the trailing direction 
and tails overlap with the orbit. It is surprising that the leading trails are not observed 
in either case. The leading trail may be a transient feature. 



• Comets 94P and 127P were classified as an "intermediate trail" by Reach et al. i.e., 
the dust more closely followed the /3 = 10~^ syndyne than the projected orbit of the 
comet (this label is consistent with our "trail" label). Neither of these comets have 
dust in our survey images and both were taken at larger true anomalies, emphasizing 
that such dust may be relatively short lived (trail survey / = 62° and 77°, SEPPCoN 
/ = 150° and -152°). 

In summary, out of the 12 comets in both surveys, three morphologies have consistent de- 
scriptions (62P, 78P, and 107P), four trails appear to be too faint for SEPPCoN (32P, 69P, 
123P, and 131P), two intermediate trails appear to be short lived and limited to smaller true 
anomalies (94P and 127P), one comet's trail may not be present at very small true anomalies 
(121P), and two trails may be obfuscated by dust tails, whereas their leading trails are possi- 
bly transient features (48P and 129P). From this comet-by-comet comparison, it appears the 
low trail detection rate of SEPPCoN can be explained by: 1) observation timing/geometry; 
2) the limited field-of-view of the IRS peak-up arrays; and, 3) the survey's sensitivity to 
dust. The last point is discussed further below. 

A portion of the difference in trail detection rates is due to the typical observing geome- 
tries of short-period comets at r/^ > 4 AU, which cause tails and trails to overlap on the sky 
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(e.g., compare 22P to 48P in Fig. [I]). However, obfuscation from brighter coma and tail dust 
does not explain the lack of trails in the 56 images of apparently bare nuclei. Instead, we 



must consider the survey differences in sensitivity. First note that the [Reach et al.| survey 
targeted comets within 3.5 AU from the Sun, whereas the SEPPCoN observations were all 
at Th > 3.4 AU. Also note that the trail survey integration times were shallow (either 42 or 
140 s on source with MIPS) compared to the SEPPCoN observations (140 to 2500 s with 
MIPS). We can make direct comparisons between the surveys by converting trail fluxes into 



optical depths. We examined the peak r values listed in Table 2 of [Reach et al.[ and find the 
trails have 0.2x10"^ < t < 9xl0~^ in the mid-IR, with a median of 1.2x10"^. For the eight 
SEPPCoN targets with trail detections, we find 0.4x10"^ < t < 11x10"^, and a median of 



1.3x10 ^. In Fig. [7| we also show a histogram of the peak optical depths from the Reach 



et al. survey. The range and median trail optical depths of the two surveys roughly agree. 



But the Reach et al. trail survey did detect most of the faintest trails, and this is likely the 



main difference between the two surveys. 



4. Discussion 



4.I. Comet activity 

Comet activity is controlled by physical and dynamical processes at the nucleus. Physical 
properties, such as composition and geology, determine if a comet is active under a given set 
of solar illumination conditions. Besides size and color temperature, little is known about the 
physical properties of our targets. The dynamical properties of a nucleus are described by its 
orbit and rotation state, and together they also determine the solar illumination conditions 
and history. The rotation states are known for only a few of our targets. In contrast, all of 
the orbital parameters are constrained well enough to make meaningful comparisons. The 
geometrical circumstances of each observation are also well constrained and can be tested. 
For example, we might test for a correlation between activity and observer-comet distance 
because dust around comets farther from the telescope is more difficult to detect (lower 
spatial resolution, possibly cooler dust temperatures). 

For this work we have tested for the presence of a correlation between recent activity 
and eight observational and orbital parameters: Vh, As {Spitzer-comet distance), (^s (phase 
angle. Sun-comet- Spitzer), f (true anomaly), q (perihelion distance), e (orbital eccentricity), 
a (orbital semi- major axis), and Ag|i5o (perihelion distance history, described below). All 
orbital parameters are measured from their osculating elements at the time of the observa- 
tion, except Ag|i5o. We integrated each comet's current orbit back 300 yr in 90 day steps 



using HORIZONS (|Giorgini et al. 1996); Ag|i5o is the difference between the minimum and 



maximum q over the past 150 yr, 

Aq = min(g) — max(g) (2) 

where a negative Ag|i5o indicates the perihelion distance was larger in the past. A relation- 
ship between any one of the above parameters and activity is not necessarily expected (e.g., 
we don't expect activity to depend on the phase angle of the observation), but we tested 
for correlations in order to perform an unbiased study while getting a useful sense of what 
a "null hypothesis" result looks like. 
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For each parameter, we compute a two-sided Kolmogorov-Smirnov (K-S) statistic, D, 
and two-tailed p-value with the null hypothesis that the active and inactive comets in our 
survey are drawn from the same distribution. The p-value (or false positive probability) 
indicates the probability that uncorrelated data sets might result yield a D greater than 
that observed. When D is small, or when the p-value is large, the null hypothesis cannot be 
rejected. At first, only those comets with comae and/or tails were considered to be active. 
Comets with trails, or without dust are considered inactive. The four comets with ambiguous 
morphologies (tail/trail) are dropped from the analysis. With these definitions, the numbers 
of active and inactive comets are 21 and 64. Our K-S test results are listed in Table HI 



In ^3.1, we argue that dust tails are associated with activity on longer timescales than 
dust comae, and that they do not necessarily indicate a currently active comet. As an 
additional check on the significance of our K-S tests, we removed tail-only detections from 
our active comet list, and compared the remaining comets to the inactive comets. This 
change reduced the number of active comets to 10. The second set of K-S tests are also 
presented in Table |4j 

The p- values in Table |4] indicate that the active and inactive comets have significantly 
different /, g, e, and a distributions. The distribution of A5 is slightly different but with 
less confidence, indicating that cov[vei- Spitzer distance does not have a strong influence on 
our results. Figure [8] presents the cumulative probability functions for /, g, a, and A^. The 
cumulative probability functions for Ag|i5o are presented in Fig. [9] We discuss each result 
below. 

4-1.1. True anomaly 

The active and inactive comets have significantly different distributions in true anomaly. 
In Fig. |8| we see that the active comets are strongly biased toward positive true anomalies 
(i.e., post-perihelion), with only 2 active comets at / < 0° (i.e., pre-perihelion) . Restricting 
our definition of active to only those comets with a coma does not diminish our conclusion. 
In addition, no comet in our survey is active between / = —180° and —125°. 

As an alternative metric of post-perihelion activity, we estimate the true anomalies at 



which activity starts and stops. Figure 10 shows a histogram of activity versus true anomaly, 
with a horizontal line plotted at exp(— 1). We chose exp(— 1) because it marks the half length 
of an exponential fall-off, however we note that our histogram does not necessarily suggest 
this particular functional relationship between activity and true anomaly. By comparing the 
line to the histogram, we estimate that activity in the average JFC turns on at / > —120° 
and turns off at / ^ 140°. 

A few scenarios may result in apparent persistent post-perihelion activity: 1) a temporary 
dust mantle insulates the ices during the pre-perihelion portion of the orbit but is lost near 
perihehon allowing persistent activity on the outbound leg of the orbit; 2) the surface layers 
retain a significant seasonal thermal wave that continues to heat the comet sub-surface when 
insolation would normally be insufficient; 3) a seasonal effect near perihelion that causes 
persistent activity well outside of perihelion; or, 4) asymmetric activity caused by the late 
onset of volatile production. For the purposes of this paper, we define a dust mantle to be 
an ice-free layer at the surface of the nucleus. 

As solar insolation warms the surface of a comet a thermal wave propagates inward. At 
some depth the thermal wave heats ices to sublimation temperatures, which drives mass loss. 
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The depth of those ices, the thermal properties of the layers above them, and the rotational 
and orbital states of the nucleus will determine when activity is initiated. Suppose that 
activity becomes so vigorous that any putative dust mantle is ejected, either partially or 
completely, from the comet. As the comet recedes from the Sun, activity would be stronger 
than the pre-perihelion orbit, but it would eventually decrease to the point at which it 
can no longer lift some grains from the surface, perhaps because of their large sizes, and 
a dust mantle will redevelop. Thus, a pre-perihelion nucleus could be insulated by a dust 
mantle, whereas the post-perihelion nucleus might not. Brin and Mendis (1979) present 
one of the first 1-D models of comet mantle formation and destruction that causes pre- and 
post-perihelion asymmetries near r^ = 1 AU. Mantle formation and destruction in nucleus 



thermal models is reviewed by Prialnik et al. (2004) 



In the second scenario, where activity is driven by seasonal heating, the dust mantle is not 
lost but retains heat from the perihelion passage. Energy deposited on diurnal timescales will 
radiate back to space on roughly the same timescale, but sub-surface temperatures can build 
up on longer timescales as a comet approaches perihelion. After the comet passes perihelion, 
the thermal wave cools to both space and the comet sub-surface. If sub-surface ices continue 
to be warmed to sublimation temperatures, despite the increasing heliocentric distance, 
activity will persist. In order to observe a pre- /post-perihelion asymmetry, the ices must be 
buried at depths closer to the annual thermal skin depth than the diurnal thermal skin depth. 
If this is true, gas production should not vary diurnally at the intermediate heliocentric 
distances in our survey, or, at least, it should be weakly correlated with rotational insolation. 
There is limited support for this scenario in flyby images of nuclei. Jet activity occurs from 
surfaces in local night at 81P/Wild 2, 9P/Tempel 1, and 103P/Hartley 2 (Sekanina et al. 



2004, Farnham et al.|[2007 , A'Hearn et al.|[20lT ), but the apparent timescales of this activity 
are much shorter than what is required for a pre-/post-perihelion asymmetry. 

The near-surface structure of comets was directly examined by the Deep Impact mission 
(A'Hearn et al. 2005). Deep Impact excavated a crater with a diameter of 50-200 m on 
Comet 9P/Tempel 1 (Richardson and Melosh 2013 Schultz et al. 2013). According to 



A'Hearn (2008), most observations of the ejecta are consistent with dust-to- volatile ratios 



of order unity, and Groussin et al. (2010) came to the same conclusion. Since the dust-to- 



volatile ratio is not ^ 1, it appears the thermal wave in this nucleus had not sublimated 
much of the water ice on the scale of the crater depth. Even greater constraints on the depth 
of thermal penetration were obtained with Deep Impact spectra of the comet's surface. 



Groussin et al. (2007) analyzed the thermal emission from the surface and showed that it 
has a low thermal inertia, nearly in instantaneous equilibrium with sunlight. They argue 
that cometary activity must originate within the first centimeters to meters of the surface. 
From the thermal inertia upper-limits and assumptions on conductivity and heat capacity. 



A'Hearn (2008) computes the diurnal skin depth to be 3 cm, and the annual skin depth to be 



90 cm. These numbers are order of magnitude consistent with theoretical models of comet 



nuclei made before the Deep Im,pact mission results (Prialnik et al. 2004), and with Spitzer 



mid-IR spectral measurements of the nucleus thermal emission at 5 AU (Lisse et al. 2005). 



In order for dust mantles to quench activity on diurnal timescales, they must be larger than 
~ 1 cm thick, and to prevent the seasonal thermal wave from driving activity, volatiles must 
be more than ~ 10 cm from the surface. 



Weissman (1987) proposed that the sudden illumination of a nucleus hemisphere caused 
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by the rapid change of season near perihehon causes cracks in the mantle, and these cracks 
create new active areas on the nucleus. The hemisphere illuminated on the in-bound leg 
of the orbit does not develop these cracks because it is gradually heated, and the surface 
can better accommodate the changes. If this hemisphere continues to be illuminated as the 
comet recedes from the Sun, the comet's activity could be greater than that of it's inbound 
leg. 



Meech et al. (2011 ) propose that the pre-/post-perihelion asymmetry of Comet 103P/Hartley 
2 is caused by the late on-set of CO2 driven activity. They find that this comet's activity 
was first initiated in 2010 at 4.3 AU by water sublimation. When the comet reached 1.4 AU 
(pre-perihelion) CO2 driven activity dominated, to the point at which water ice can be driven 
off the surface, as if it were dust. The CO2 driven activity would persist beyond the comet's 



aphelion ( [Meech et al.||201lD . 

Whatever the cause, the effect is clear in the SEPPCoN. Because the survey took a 
consistent approach to its observations, it can be used to test these and other inhibitors and 
drivers of activity in a statistical manner. 

4- 1-2. Perihelion distance 

The inactive comets have smaller perihelion distances than the active comets. The median 
q is 2.0 AU for the inactive comets and 3.0 AU for the active comets. This difference may 
be a selection effect. Cometary activity is driven by insolation, and, in general, comets are 
most likely to be active when they are near perihelion. Our survey targets were observed at 
heliocentric distances ranging from 3 to 7 AU. However, a comet with a perihelion distance 
between 3 and 7 AU is a priori more likely to be active in our survey, since its activity at 
such distances has already been established (otherwise it would not have been designated a 
comet, and not observed by SEPPCoN). Therefore, the greater rate of activity for comets 
with larger perihelion distances appears to partially be an artifact of the survey. 

To further demonstrate the larger activity rates of the large perihelia comets, we plot 
the cumulative frequency of activity in our survey, sorted by perihelion distance, in Fig. [TT| 
That the frequency of activity almost continuously increases for q > 2 AU shows the higher 
activity rate of this population. 



Figure 12 presents an alternative view on these large perihelia comets, and shows that 
on average they produce more dust than the smaller perihelia comets. Here, we plot efp 
versus perihelion distance. We define the parameter efp as the product of effective emissivity 
(e), filling factor (/), and projected aperture radius (p, in units of cm). This parameter is 
intended to be analogous to the Afp parameter commonly computed for scattered light 



observations (A'Hearn et al. 1984). Further discussion of efp can be found in Appendix 
\N\ We based our efp values on the mean MIPS 24-iim, or the IRS 22-iim center aperture 
photometry (Table pi). Comets with small efp values, i.e., low-activity comets, are missing 
from the large-perihelion population. We suggest that this bias is due to the discovery 
circumstances of large-perihelia comets, and that low-activity comets with large-perihelia 
are as yet missing from the discovered JFC population. 

Although we cannot make definite conclusions on the activity of comets with large- 
perihelia, we can consider the activity of small-perihelia comets, whose discovery rates are 
likely more complete and less biased. Of those comets with q < 2.5 AU, 7 out of 59 comets 
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Figure 11: The cumulative frequency of activity in our survey versus perihelion distance. 
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Figure 12: (Left) Comet activity, quantified with the parameter efp, versus perihehon distance. (Right) The 
same as on the left, but efp has been scaled by {rfi/3 AU)^, so that comets can be more directly compared 
to each other. Note the paucity of low-e/p comets with large perihelion distances, indicating that our survey 
and the discovered JFC population is sensitivity limited. 
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(12%) are active, and none of the 30 comets with q < 1.8 AU are active (Fig.[TT|p] Therefore 
any short-period comet with a small perihelion distance is likely inactive at 3-7 AU. This 
finding may be explained if relatively more volatile ices are lost for comets that approach 
the Sun more closely during their perihelion passages, but this hypothesis seems inconsistent 
with the results of A'Hearn et al. (2012). They studied the H2O, CO2, and CO content 



of comets in the literature, and found no correlation with orbital parameters, except for 
a decreasing C02-to-H20 mixing ratio for comets with decreasing perihelion distance, for 
g < 1 AU. They rightly dismiss this correlation as due to small number statistics (only 4 
comets have a measured CO2 abundance for g < 1 AU) and for the fact that it contradicts 
their result that CO/H2O is uncorrelated with q. Taking our results together with the 
discussion of A'Hearn et al. (2012), we suggest that comets with small perihelion distances 



have a thicker insulating surface than comets with larger perihelion distances. 
4.. 1.3. Perihelion distance history 



Licandro et al. (2000) surveyed 18 comets in order to estimate the effective sizes of 
their nuclei. Seven of their targets were active: one at r^ = 3.04 AU, and the rest at 
Th > 4: AU. They examined their orbital histories and found that all but one of their comets 
recently perturbed to orbits with smaller perihelion distances in the past 150 yr (Ag|i5o < 



—1 AU) were active. Licandro et al. hypothesize that the perturbations to smaller perihelion 
distances caused a marginally stable mantle to be destroyed, exposing fresh ices. The absence 
of a mantle allows activity to persist to larger than usual perihelion distances. When we 
examined the pre-/post-perihelion asymmetry in SEPPCoN, we also introduced the concept 
of a temporary mantle, only in our case the mantle is partially or completely lost each 
perihelion passage. We now examine the orbital histories of SEPPCoN targets to further 
investigate the potential existence of dust mantles. Note that a small Ag is not the same 
as small q, i.e., a comet may be perturbed from the Centaur region (a > 5.2 AU) into the 
Jupiter-family comet region (a < 5.2 AU), yielding Ag < 0, yet still have a large perihelion 
distance in our survey. 

Recognizing that the same Ag is more significant at small perihelion distances than 
at large perihelion distances, we have also computed the maximum sub-solar-temperature 
difference, AT. The sub-solar temperature of a spherical comet nucleus is 



T 



A)*Fq 



nv^f^sB 



1/4 



(3) 



where A is the Bond albedo, Fq = 1365 W m ~^ is the solar constant at 1 AU, r^ is the 
heliocentric distance in AU, 77 = 1.03 is the IR-beaming parameter, e = 0.95 is the IR 
emissivity of the surface, and asB is the Stefan-Boltzmann constant in W m~^ K"^ (Harris 



1998, Fernandez et al. 2011[ ). For our purposes, we compute AT\x using the minimum and 
maximum g values from our Agl^, integrations; AT > signifies the temperature was cooler 
in the past. 



^Because the strong bias in activity to post-perihelion epochs, we verified that the 30 comets with q < 
1.8 AU have the same true anomaly distribution as comets with q > 1.8 AU. A K-S test yields a 1% 
probability that the two populations are different. 
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Figure [9] shows that active comets in SEPPCoN appear more hkely to have Ag|i5o < and 
AT|i5o > than the inactive comets. In addition, no active comet has Ag|i5o > 0.5 AU, yet 
~ 10% of the inactive comets do have such large values. The K-S probabilities for Ag|i5o and 
AT|i5o are 28% and 10%, indicating that the difference between active and inactive comets 
may be significant. But, given the striking differences in Fig. [9} this problem deserves future 
study because it may lead to an understanding of the timescales of mantle formation. 

The difference in perihelion distance histories between the active and inactive comets 
may be related to their discovery circumstances. When an unknown comet is perturbed 
to a smaller perihelion distance, it should become brighter to observers on the Earth, and 
therefore more likely to be discovered. To investigate this possible effect, we searched the 
NASA JPL Small-Body Database Browser and the NASA Planetary Data System Small 
Bodies Node for each comet's year of discovery. We plot this year versus the year of the 



minimum Ag|i5o in Fig. 13 Data points near a slope of 1 indicate a comet that was discovered 
at a 150-year low in its perihelion distance history. Fifteen comets are found within 20 years 
of this line, indicating that reductions in perihelion distance is a factor in JFC discovery 
dates. In the plot we indicate which comets were active in our survey. Those comets 
discovered soon after being scattered to smaller perihelion distances appear to have a higher 
activity rate in our survey, 8 out of 22 (36 ± 13%) versus 15 out of 70 (21 ± 6%), but the 
low significance precludes us from making a firm conclusion. 

Most of our survey targets were discovered in the past 20 years, but the year of minimum 
Ag|i5o for these comets is uniformly distributed over the past 150 years. This distribu- 
tion may be the result of uncertain orbital histories for these comets, rather than a lack 
of correlation between the two parameters. However, it is also possible that the increasing 
sophistication and sensitivity of Solar System surveys in the past 20 years has led to JFC dis- 
coveries that are independent of their perihelion distance history. Both of these possibilities 
could be a cause for the low significance of the K-S tests based on Ag. 

4-1 .4- Semi-major axis and Spitzer-com,et distance 

On average, the active comets have larger semi-major axes than the inactive comets. 
Since perihelion distance and semi-major axis are correlated in the Jupiter-family comet 
population, this difference may be an observational bias related to the perihelion distance 
bias described above. A similar argument may be made for the slight differences in A5 for 
the active and inactive comets. Nuclei with larger q are more likely to be observed at larger 
A5. 

4.1.5. Other investigations 

Comet lightcurve asymmetries have been recognized for many years. Asymmetries near 
perihelion are not directly relevant to our survey since most of our targets are observed at 
true anomalies greater than 90° from perihelion. The secular light curves of 27 comets were 



studied by Ferrin (2010). Considering only the Jupiter- family comets, and excluding the one- 
time active Comet 107P/Wilson-Harrington (leaving 17 comets), the ratios of the heliocentric 
distance at which activity ceases (^o//) to the distance at which activity initiates (r^^) ranges 
from 0.7 to 2.0, with a mean Toff/ron = 1-20 ± 0.38. Although this result cannot be directly 
compared to our survey results, they qualitatively agree in that it appears persistent post- 
perihelion activity is common to individual Jupiter-family comets. 
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Figure 13: Year of comet discovery iodise) versus the year of the comet's 150-year minimum in pcriheUon 
distance history (TaijIiso)- Comets between the dotted and sohd trend hues are those that were discovered 
near a 150-year minimum in pcrihehon distance: 8 out of the 22 comets in this region are active in our 
survey. 
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The activity of short-period comets has also been examined by Mazzotta Epifani et aL 



(2009). They searched the hterature for optical and infrared observations of short-period 
comets (Jupiter-family, Encke type, and Halley type) at r^j > 3 AU. In their compilation 
of 90 comets they find that 10-20% are active beyond 4.5 AU. Furthermore, they also find 
a perihelion asymmetry with activity detected in 22% of the pre-perihelion observations, 
and 59% of the post-perihelion observations. In our survey data {rt > 3.4) we find slightly 
different values: 24±5% (21 comets) of our sample was active, with activity in 5±4% (2/39) 
of the pre-perihelion observations, and 38 ± 9% (19/50) in the post-perihelion observations. 



With only a few comets in each bin, and understanding that the Mazzotta Epifani et al. 



study is based on a literature search, and not on a consistent data set, we consider the 



SEPPCoN results to be compatible with the results of [Mazzotta Epifani et aLj In addition, 
they also conclude that recent perihelion distance changes are not a strong predictor of 
distant activity, agreeing with our results above. 

4-2. Dust color temperature 

The effective temperatures of comet comae at 5-20 pm typically range from ?^ — 30% 



warmer than an isothermal blackbody sphere at the same heliocentric distance (Gehrz and 
Ni^[l992| psse et al.|[l998| [Sitko et"aL][2004[ [Woodward et aL||20lI| . In Table |3j we find 



a similar result, even though our observations are obtained at larger heliocentric distances 
and longer wavelengths than are typically measured with ground-based observatories. 

We test for correlations between Tc/Tbb for the 11 IRS center aperture measurements and 
the parameters rh, q, -F22 (22- or 24-iim flux), and efp. The results are presented in Table [5J 
For each parameter, we calculate the linear correlation coefficient and apply a Spearman "p" 
test. The Spearman test is a more robust estimate of correlation than the linear correlation 
coefficient. It is based on the relative rankings of each parameter being tested, and does not 
test for a specific functional form of the correlation. After computing the Spearman p, we 
derive Z, the number of standard deviations by which the Spearman statistic deviates from 
the null hypothesis expectation value for uncorrelated data. To account for measurement 
uncertainties in the correlation tests, we tested 10'^ data sets generated with a Monte Carlo 
technique based on our measurements and associated uncertainties. No significant correlation 
is found, and the correlations in our Monte Carlo runs never exceeded a 3a significance. We 
list our results in Table M 

5. Conclusions 

We presented an analysis of the dust and activity of comets observed by the Spitzer Space 
Telescope as part of the Survey of the Ensemble Physical Properties of Cometary Nuclei 
(SEPPCoN). We have shown that SEPPCoN can be used to study the dust and activity 
of Jupiter-family comets at 3-7 AU, although some care must be taken when interpreting 
the results as the dust detections are limited by the image sensitivities, and that comets 
with large perihelion distances are a priori more active. We detected dust around 33 of 89 
(37 ± 6%) survey targets, and 21 targets (24 ± 5%) have comae or tails that suggest recent 
cometary activity. Since faint dust (r < 0.5x10^^) may remain undetected in our images, we 
conclude our detections are a lower limit and that at least ~ 24% of Jupiter-family comets 
are active at 3-7 AU from the Sun. We also studied the crude spectral properties of the 



dust. The 16- to 22-]am color temperature of the dust is 7.4 ± 0.6% warmer on average than 
an isothermal blackbody sphere in LTE. 

We introduce the quantity efp, intended to be a thermal emission counterpart to the 
often reported Afp for observations of light scattered by comets. The efp versus perihelion 
distance distribution of our survey sample shows that low-activity comets with large perihe- 
lion distances are missing from the survey, and therefore are likely missing from the known 
Jupiter-family comet population. 

We compared the frequency of activity of survey targets to their orbital and observational 
parameters. All 30 comets with g < 1.8 AU appear inactive in our survey. Whatever drives 
their activities is effectively quenched at larger distances from the Sun. 

Comets that have been perturbed to smaller perihelion distances in the past 150 years 
seem more likely to be active in our survey (cf. [Licandro et al. 2000). Kolmogorov-Smirnov 



tests indicate a tentative result, so a better statistical characterization, perhaps with more 
detailed orbital histories, is warranted. 

The activity of Jupiter-family comets at 3-7 AU from the Sun is significantly biased to 
post-perihelion epochs. Of the 21 comets that appear to be recently active, 19 were observed 
post-perihelion. We tentatively estimate the true anomaly at which activity is initiated to be 
/ > —120°, but small number statistics make this limit uncertain. With better confidence, 
we find that for many comets activity has ceased by / ~ 140°. Similar findings on the 



persistent post-perihelion activity were found by Mazzotta Epifani et al. (2009), and Ferrin 



(2010). No comet in our survey is active between / = —180° and —125°. 

We discussed several causes for a pre-/post-perihelion activity asymmetry. Whatever the 
cause, the effect is clear in our survey. Because SEPPCoN took a consistent approach to 
its observations, it can be used to test the drivers and inhibitors of activity in a statistical 
manner. The Rosetta spacecraft, which will orbit Comet 67P/Churyumov-Gerasimenko 



(Glassmeier et al. 2007) as it approaches perihelion from 4 AU, should also provide new 



perspectives on this problem of the evolutionary history of comet surfaces. 
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Tabic 1: SEPPCoN target list, observational geometry, and orbital parameters: Obs. date — Spitzer observation date, 
for MIPS observations, only one of the two epochs are listed; rj, — heliocentric distance; Ag — Spitzer comet distance; 
05 - Snn-comet- Spitzer angle; / — true anomaly, 0° at perihelion; T — Tp - time from perihelion, < when / < 0; a 
— semi-major axis; e - eccentricity; q - perihelion distance; Ag|i5o - change in q over the past 150 yr. 



Name 



Obs. date 
(UT) 



rh 
(AU) 



As 
(AU) 



n 



f 

(°) 



T -Tp a 

(days) (AU) 



(AU) 


A<j|l50 

(AU) 


1.35 


-0.23 


1.25 


0.50 


1.56 


0.56 


2.72 


1.25 


0.97 


0.21 


1.47 


0.43 


1.58 


-0.57 


3.42 


1.57 


1.83 


-0.42 


2.17 


0.78 


1.57 


-0.19 


2.79 


-1.76 


2.31 


-0.35 


1.92 


-0.15 


1.68 


0.19 


2.53 


0.12 


1.72 


0.46 


1.48 


-0.84 


1.76 


-0.29 


1.94 


0.42 


3.56 


-2.14 


2.31 


-0.90 


2.01 


-1.01 


1.23 


-0.22 


2.28 


-0.30 


1.71 


0.07 


2.24 


-0.43 


2.35 


-0.31 


0.99 


0.05 


2.13 


0.36 


2.00 


-0.21 


3.04 


-2.24 


2.65 


-1.10 


2.13 
1.47 


-0.23 
-0.45 


2.18 


-0.21 


2.81 


-3.06 


2.10 


-0.21 


2.42 


-1.24 


1.92 


0.28 


1.89 


0.21 


1.71 


0.14 


3.40 


-0.83 


2.54 


0.65 


1.44 


-0.11 


1.39 


0.15 


1.70 


-1.36 


2.65 


-0.70 


3.11 


-1.28 


3.65 


-0.14 


2.08 


-0.68 


1.23 


0.14 


1.92 


0.27 


1.42 


1.17 


0.61 


-0.05 


1.73 


-0.84 


2.24 


-1.09 


4.21 


0.09 



6P/d'Arrest 2007-02-13 22:48 

7P/Pons-Winnecke 2007-04-19 14:39 

llP/Tempel-Swift- LINEAR 2006-09-01 15:12 

14P/Wolf 2007-03-20 17:52 

15P/Finlay 2007-02-28 03:06 

16P/Brooks 2 2007-04-10 09:32 

22P/Kopff 2007-04-19 14:10 

31P/Schwassmann-Wachmann 2 2006-11-09 10:30 

32P/Comas Sola 2006-08-01 03:59 

33P/Daniel 2006-11-09 17:14 

37P/Forbes 2007-03-12 00:56 

47P/Ashbrook-Jackson 2007-03-21 10:56 

48P/ Johnson 2006-11-11 18:43 

50P/Arend 2006-10-23 18:43 

51P-A/Harrington 2007-04-06 23:24 

56P/Slaughter-Burnham 2006-12-23 13:48 

57P-A/du Toit-Neujmin-Delaporte . . 2007-06-15 08:34 

62P/Tsuchinshan 1 2006-10-03 02:00 

68P/Klemola 2007-02-12 09:42 

69P/Taylor 2006-07-29 05:19 

74P/Smirnova-Chernykh 2007-01-28 17:01 

77P/Longmore 2007-02-10 07:38 

78P/Gehrels 2 2007-03-10 02:00 

79P/du Toit-Hartley 2006-09-17 23:10 

89P/Russell 2 2007-05-01 09:31 

93P/Lovas 1 2006-10-03 23:46 

94P/Russcll 4 2006-11-14 19:47 

lOlP-A/Chernykh 2007-04-22 13:51 

107P/Wilson-Harrington 2007-02-12 07:43 

113P/Spitaler 2006-10-22 15:58 

118P/Shocmaker-Levy 4 2006-09-17 11:45 

119P/Parker-Hartley 2007-02-12 08:01 

121P/Shocmaker-Holt 2 2006-08-01 02:40 

123P/West-Hartley 2006-10-21 23:31 

124P/Mrkos 2006-09-16 19:19 

127P/Holt-01mstead 2007-03-08 16:37 

129P/Shocmaker-Levy 3 2007-04-19 00:39 

130P/McNaught-Hughes 2006-11-14 10:57 

131P/Mueller 2 2007-02-12 08:54 

132P/Helin-Roman-Alu 2 2007-06-07 19:26 

137P/Shocmaker-Levy 2 2007-03-16 13:11 

138P/Shoemaker-Levy 7 2007-01-22 19:49 

139P/Vaisala-Oterma 2006-11-02 07:11 

143P/Kowal-Mrkos 2007-03-09 16:14 

144P/Kushida 2007-07-09 22:38 

146P/Shoemaker-LINEAR 2006-08-06 20:20 

148P/Anderson-LINEAR 2006-11-02 10:30 

149P/Mucller 4 2007-02-09 03:36 

152P/Helin-Lawrence 2006-09-17 20:19 

159P/LONEOS 2007-02-13 08:58 

160P/LINEAR 2006-12-19 16:49 

162P/Siding Spring 2007-03-17 23:43 

163P/NEAT 2006-08-06 17:59 

168P/Hergenrother 2007-05-23 02:26 

169P/NEAT 2007-03-01 18:40 

171P/Spahr 2007-03-16 13:00 

172P/Yeung 2006-10-18 20:36 

173P/Mueller 5 2006-10-24 09:18 

Continued on next page 



4.39 3.85 12 -145 -548 3.49 0.61 
4.34 4.03 13 -146 -526 3.43 0.63 

4.40 4.04 13 -146 -611 3.42 0.54 
4.52 4.53 13 -120 -709 4.25 0.36 
4.40 4.32 13 -149 -480 3.48 0.72 
3.40 3.37 17 -125 -368 3.36 0.56 
4.87 4.38 11 -156 -767 3.46 0.54 
5.02 4.59 11 -167 -1423 4.23 0.19 
4.14 3.76 14 122 486 4.27 0.57 

4.40 4.17 13 -127 -619 4.03 0.46 

4.31 4.19 13 144 587 3.43 0.54 

4.32 4.32 13 -117 -682 4.12 0.32 

4.41 4.32 13 141 761 3.65 0.37 

3.48 3.31 17 -107 -373 4.09 0.53 
3.77 3.35 15 -125 -438 3.70 0.55 
5.08 5.04 12 120 708 5.10 0.50 

4.11 3.92 14 -138 -560 3.45 0.50 
4.72 4.66 12 150 664 3.52 0.58 
5.47 5.28 10 -138 -709 4.89 0.64 
4.25 3.66 12 135 605 3.65 0.47 

4.42 4.01 12 -121 -914 4.17 0.15 
4.59 4.18 12 -152 -878 3.60 0.36 

4.98 4.46 10 153 865 3.74 0.46 
4.37 4.08 13 -158 -618 3.03 0.59 
4.76 4.18 11 -146 -839 3.80 0.40 
4.01 3.56 14 -121 -439 4.39 0.61 
4.79 4.20 11 178 1174 3.51 0.36 

4.34 3.76 12 103 483 5.78 0.59 
4.08 3.52 13 166 582 2.64 0.62 
3.86 3.55 15 -120 -518 3.69 0.42 
4.95 4.58 12 178 1159 3.47 0.42 
4.21 3.67 12 103 629 4.29 0.29 

4.35 3.97 13 123 696 4.02 0.34 
5.34 4.91 10 161 1049 3.86 0.45 
4.23 3.76 13 -149 -588 3.21 0.54 
4.59 4.43 13 -164 -957 3.44 0.37 
4.01 3.67 14 119 679 3.76 0.25 
4.45 4.20 13 146 752 3.54 0.41 
4.42 3.96 12 140 787 3.68 0.34 

3.99 3.80 15 120 478 4.09 0.53 
5.47 5.45 11 -142 -788 4.48 0.58 
4.21 4.21 14 136 552 3.63 0.53 
4.10 3.57 13 -82 -535 4.51 0.25 
4.99 4.74 11 -133 -825 4.30 0.41 
4.61 4.16 12 -142 -567 3.86 0.63 

5.08 4.86 12 -147 -649 4.03 0.66 
4.31 3.89 13 -137 -567 3.68 0.54 
5.55 5.46 10 -150 -1106 4.33 0.39 
5.70 5.58 10 158 1364 4.50 0.31 

6.12 5.74 9 118 1078 5.88 0.38 
4.99 4.83 12 144 798 3.98 0.48 
4.82 4.27 11 173 858 3.05 0.60 

4.09 3.95 14 130 552 3.67 0.48 

4.49 4.09 12 144 567 3.63 0.61 
4.29 4.01 13 168 530 2.60 0.77 
4.16 4.09 14 137 559 3.53 0.51 
4.25 4.06 14 -141 -725 3.51 0.36 
4.82 4.32 11 -67 -572 5.71 0.26 
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Table 1: continued 



Name 



Obs. date 
(UT) 



(AU) 



As 
(AU) 



<Ps 

n 



f 

(°) 



T -Tp a 

(days) (AU) 



(AU) 



A<j|l50 

(AU) 



197P/LINEAR 2006-11-04 17:12 

213P/2005 R2 (Van Ness) 2006-11-13 12:57 

215P/2002 OS (NEAT) 2007-01-08 17:45 

216P/2001 CVg (LINEAR) 2007-01-01 12:48 

219P/2002 LZii (LINEAR) 2007-06-07 04:31 

221P/2002 JNi6 (LINEAR) 2007-04-10 02:54 

223P/2002 SI (Skiff) 2007-04-06 00:18 

228P/2001 YXi27 (LINEAR) 2006-10-22 11:26 

246P/2004 F3 (NEAT) 2006-12-17 23:20 

256P/2003 HT15 (LINEAR) 2006-12-06 12:30 

260P/2005 K3 (McNaught) 2007-01-01 07:44 

C/2005 W2 (Christenscn) 2007-01-01 04:24 

P/2001 R6 (LINEAR-Skiff) 2007-04-05 14:06 

P/2003 03 (LINEAR) 2007-07-09 12:00 

P/2004 Al (LONEOS) 2007-03-05 01:21 

P/2004 DO29 (Spacewatch-LINEAR) 2006-09-18 04:24 

P/2004 H2 (Larsen) 2007-01-01 20:40 

P/2004 V3 (Siding Spring) 2006-09-18 05:33 

P/2004 V5-A (LINEAR-Hill) 2007-03-10 01:13 

P/2004 VRg (LONEOS) 2006-08-07 18:51 

P/2005 GFg (LONEOS) 2006-11-12 23:27 

P/2005 JD108 (Catalina-NEAT) 2007-01-01 21:20 

P/2005 JQ5 (Catalina) 2006-11-21 13:58 

P/2005 L4 (Christensen) 2007-03-11 23:02 

P/2005 Q4 (LINEAR) 2007-02-28 09:05 

P/2005 Rl (NEAT) 2006-12-19 12:41 

P/2005 S3 (Read) 2007-05-23 02:48 

P/2005 T5 (Broughton) 2006-10-09 13:43 

P/2005 W3 (Kowalski) 2006-12-08 04:36 

P/2005 XA54 (LONEOS-Hill) 2007-06-04 23:12 

P/2005 Y2 (McNaught) 2007-01-20 05:02 



4.21 


4.02 


14 


-159 


-561 


2.87 


0.63 


1.06 


-0.15 


4.05 


3.71 


14 


137 


641 


3.43 


0.38 


2.13 


-0.22 


4.77 


4.24 


11 


-161 


-1245 


4.03 


0.20 


3.22 


-1.01 


4.41 


3.95 


12 


-131 


-648 


3.89 


0.44 


2.16 


-0.43 


4.78 


4.22 


11 


-160 


-1004 


3.65 


0.35 


2.37 


0.12 


4.40 


3.83 


12 


-145 


-656 


3.48 


0.49 


1.79 


-0.31 


5.70 


5.15 


9 


-163 


-1229 


4.14 


0.42 


2.41 


-0.74 


4.84 


4.44 


12 


162 


1330 


4.16 


0.18 


3.43 


-0.78 


4.28 


3.71 


12 


119 


712 


4.02 


0.29 


2.87 


-1.01 


6.21 


5.66 


-8 


158 


1328 


4.61 


0.42 


2.68 


-0.22 


4.20 


3.94 


14 


136 


508 


3.69 


0.59 


1.51 


0.20 


4.07 


3.76 


14 


53 


280 


18.93 


0.82 


3.33 


-0.02 


5.74 


5.18 


9 


-158 


-1086 


4.15 


0.48 


2.17 


-0.29 


4.28 


4.10 


14 


-153 


-571 


3.10 


0.60 


1.25 


-0.60 


6.52 


6.42 


9 


71 


910 


7.96 


0.31 


5.48 


-4.38 


5.55 


5.21 


10 


82 


712 


7.36 


0.44 


4.09 


-0.14 


5.46 


5.07 


10 


140 


965 


4.51 


0.42 


2.62 


0.23 


5.36 


4.90 


10 


82 


676 


7.12 


0.45 


3.94 


-0.20 


5.81 


5.57 


10 


77 


737 


7.98 


0.45 


4.42 


-0.03 


3.44 


3.35 


17 


85 


339 


4.86 


0.51 


2.38 


-0.19 


4.16 


3.79 


14 


87 


452 


5.85 


0.52 


2.83 


0.24 


4.77 


4.32 


11 


65 


509 


6.45 


0.38 


4.03 


-0.90 


4.02 


4.00 


15 


160 


482 


2.69 


0.69 


0.83 


-0.06 


4.15 


4.07 


14 


116 


564 


4.12 


0.43 


2.37 


0.76 


4.41 


3.93 


12 


127 


518 


4.46 


0.61 


1.75 


-0.12 


4.12 


3.72 


14 


108 


-276 


5.50 


0.63 


2.05 


0.16 


4.12 


3.88 


14 


93 


500 


4.90 


0.42 


2.84 


-0.15 


3.99 


3.71 


15 


62 


340 


7.26 


0.55 


3.25 


-0.09 


4.36 


4.35 


13 


84 


472 


6.40 


0.53 


3.01 


-0.24 


4.42 


3.86 


12 


116 


454 


6.13 


0.71 


1.77 


0.11 


5.41 


5.07 


10 


101 


753 


6.30 


0.47 


3.36 


0.14 
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00 



Table 2: Dust morphology and photometry. Table columns: Instrument - the name of the instrument used; Morph. ~ Morphology "N" for no dust 
detected, "C" for coma, "Ta" for tail, "Tr" for trail, a "?" indicates a tentative detection of dust; p ~ the radius of the circular aperture centered 
on the comet nucleus; For MIPS, Fi and F2 are the flux densities of the two 23.7-]jm observations; For IRS, Fi and F2 are the flux desnities of the 
15.8-]jm and 22.3-]jm peak up arrays; d - The distance of the tail/trail aperture from the center of the comet nucleus; A - The area of the aperture 
used to measure the tail/trail flux. 



Center aperture 



Tail or trail aperture 



Comet 



Instru- 
ment 



Morph. 



P 

(Px) 



(mJy) 



F2 (72 

(mJy) 



d 

(px) 



A 
(px) 



^1 Ti 

(mJy) 



F2 02 

(mJy) 



6P/d'Arrest IRS Tr? 6 

16P/Brooks 2 MIPS Ta/Tr 6 0.528 0.062 

22P/Kopff IRS Tr 6 0.726 0.075 1.16 

32P/Comas Sola IRS C+Ta 6 16.21 0.27 33.10 

37P/Forbcs IRS Ta/Tr 6 0.625 

48P/Johnson IRS Ta 6 0.61 0.16 0.69 

50P/Arend IRS Ta? 

56P/Slaughter-Burnham IRS Ta/Tr 3 

62P/Tsuchinshan 1 MIPS Tr 4 

69P/Taylor IRS C 4 

74P/Smirnova-Chernykh IRS C+Tr 6 

78P/Gehrels 2 IRS Ta/Tr 

lOlP-A/Chcrnykh IRS C 6 

118P/Shocmaker-Lcvy 4 IRS Ta 6 

119P/Parker-Hartley IRS Ta 6 

121P/Shocmaker-Hoh 2 IRS Tr 

129P/Shoemaker-Levy 3 IRS Ta 6 

144P/Kushida MIPS Tr 6 

152P/Helin-Lawrence IRS Ta 6 

159P/LONEOS IRS Ta 6 

171P/Spahr IRS Tr 

173P/Mucller 5 IRS Ta 6 6.75 0.26 

213P/2005 R2 (Van Ness) IRS C+Ta 6 14.96 0.39 

219P/2002 LZn (LINEAR) .... MIPS Tr 4 

246P/2004 F3 (NEAT) IRS C+Ta 6 

260P/2005 K3 (McNaught) .... MIPS Tr 

C/2005 W2 (Christensen) MIPS Ta 6 

P/2004 Al (LONEOS) MIPS C+Ta 6 6.65 0.14 

P/2004 H2 (Larsen) MIPS C 4 0.415 0.062 

P/2004 V5-A (LINEAR-Hill) . . IRS Ta 6 0.640 0.048 

P/2004 VRg (LONEOS) IRS Ta 6 

P/2005 JD108 (Catalina-NEAT) MIPS C+Ta 6 10.20 0.14 

P/2005 T5 (Broughton) MIPS Tr 4 

P/2005 W3 (Kowalski) MIPS Ta 5 2.30 0.14 

P/2005 Y2 (McNaught) MIPS C+Ta 6 13.94 0.22 



0.61 0.16 

0.177 0.034 

6.52 0.19 

1.136 0.060 

0.248 0.074 

1.81 0.10 



0.315 0.036 



0.242 
0.477 
13.30 

3.20 

4.16 

4.65 

0.427 
0.539 

16.96 
32.35 



25.57 0.48 59.2 



1.41 
6.10 

1.463 
11.6 

10.53 
0.44 
2.27 

15.79 



0.17 

0.55 

0.051 

0.21 



0.050 

0.031 

0.19 

0.13 

0.18 

0.23 

0.072 
0.060 

0.27 
0.74 

1.1 

0.14 
0.11 

0.087 
1.1 
0.14 
0.12 
0.15 
0.23 



5 

11 

5 

7 

12 

15 
12 

14 
13 

6 

12 
13 

7 
27 
12 
15 
15 
15 
14 
37 
11 
38 
24 
21 

14 
12 
23 

22 
12 
17 



16.9 
91.8 
138.9 
58.3 
47.9 

77.2 
72.0 

104.2 
59.8 

64.9 

88.0 

148.4 

43.6 

142.2 

57.4 

57.6 

50.2 

111.9 

224.7 

217.9 

327.5 

402.2 

163.9 

217.9 

79.3 
125.5 
212.5 
129.1 

81.9 
173.1 



0.180 
1.109 

8.21 
0.205 

0.63 

0.216 
0.377 

1.55 
0.320 



1.914 
2.33 



0.321 

0.259 

0.285 

4.65 

21.82 

1.01 

36.33 

2.33 

2.04 

6.55 

0.350 
14.87 
9.68 
1.90 
1.85 
6.79 



0.023 
0.067 

0.30 
0.028 

0.10 

0.027 
0.040 

0.18 
0.035 



0.088 
0.13 



0.038 
0.025 
0.030 
0.26 
0.57 
0.11 
0.88 
0.12 
0.18 
0.20 

0.040 
0.79 
0.20 
0.25 
0.14 
0.27 



2.53 
19.15 
0.536 

0.99 



3.28 
0.804 

0.405 

4.27 

4.98 

3.05 

0.61 

0.818 

0.891 

0.293 

13.03 

45.6 

0.98 

87.6 

2.68 
5.98 



0.' 

28.8 
10.13 

1.58 



0.16 

0.62 

0.036 

0.13 



0.315 0.058 



0.18 
0.084 

0.085 

0.14 

0.23 

0.14 

0.18 

0.050 

0.042 

0.038 

0.27 

1.1 

0.10 

2.0 

0.17 
0.16 

0.072 

1.2 

0.20 

0.20 



7.37 0.29 



Table 3: Color-temperature excess (Tc/Tbb) for IRS peak up observations with photometry in both filters 
(Table [2]). Only temperatures with uncertainties less than 40 K are listed. Table columns: Tbb, the 
temperature of an isothermal blackbody sphere in local thermodynamic equilibrium at the same r^', Center, 
the color temperature and error, normalized by Tbbi for the center aperture of Table [2] Tail or Trail, the 
same as Center, but for the tail or trail aperture. 



Comet 



Th 



T, 



BB 



Center 



(AU) (K) TJT, 



BB 



a 



22P 4.87 126 

32P 4.14 137 

37P 4.31 134 

48P 4.41 132 

74P 4.42 132 

78P 4.98 125 

lOlP 4.34 133 

119P 4.21 135 

121P 4.35 133 

152P 5.70 116 

159P 6.12 112 

173P 4.82 127 

213P 4.05 138 

246P 4.28 134 

P/2004 V5-A 5.81 115 

P/2004 VRg . 3.44 150 



1.43 
1.12 



1.16 

0.96 
1.06 



1.52 
1.08 
1.07 
1.06 
1.24 



0.23 
0.02 



0.02 

0.03 
0.04 



0.17 
0.02 
0.02 
0.01 
0.07 



Tail or Trail 



TjTj 



BB 



1.14 

1.04 
1.00 
1.37 
1.13 
1.09 

1.07 

1.12 
1.16 
1.04 
1.02 
1.09 
1.04 
1.11 
1.05 



a 



0.05 
0.03 
0.08 
0.24 
0.09 
0.09 

0.03 
0.05 
0.08 
0.05 
0.03 
0.02 
0.02 
0.07 
0.04 



Table 4: Kolmogorov-Smirnov statistic, D, and p- value used to determine if active comets and inactive comets 
are drawn from different distributions, based on 7 parameters: r^, heliocentric distance; A5, Spitzer-coraet 
distance; ^5, phase angle; /, true anomaly; a, orbital semi-major axis; e, orbital eccentricity; g, perihelion 
distance; Agjiso, change in q over the past 150 years. Active comets are those 21 targets that were observed 
to have a coma and/or tail. As an additional test, we removed comets with tail-only morphologies from the 
active list, leaving 10 comets with coma detections. 





Active vs. Inactive 


Coma vs. Inactive 


Parameter 


D 


P 


D 


P 


Th 


0.18 


0.67 


0.19 


0.88 


A5 


0.29 


0.12 


0.36 


0.17 


05 


0.20 


0.52 


0.33 


0.24 


/ 


0.51 


0.00024 


0.48 


0.025 


a 


0.45 


0.0021 


0.43 


0.053 


e 


0.34 


0.039 


0.38 


0.13 


Q 


0.50 


0.00045 


0.44 


0.051 


Aq 


150 


0.24 


0.28 


0.31 


0.33 


AT 


150 


0.30 


0.10 


0.26 


0.52 
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Table 5: efp values for all center aperture photometry in Table [2] e/pj^g and ^f P22 ^^^^ computed from 
IRS peak-up photometry, ef P2i is computed from MIPS photometry. For comets with two epochs of MIPS 
photometry, e/p is computed from the average. For all comets, e is computed assuming a grain temperature 
10% warmer than a blackbody in LTE. 



Comet 




e/Pi6 


0"16 


e/P22 


0"22 


e/P24 ^24 






(cm) 


(cm) 


(cm) 




16P ... 












1.4 0.2 


22P ... 




20.1 


2.1 


13.2 


2.0 








32P ... 




230.3 


3.8 


223.6 


3.7 








37P ... 








5.1 


0.4 








48P ... 




12.1 


3.2 


6.2 


1.9 








62P ... 














.3 C 


.3 


69P ... 








5.0 


0.3 








74P ... 




120.8 


3.6 


110.7 


1.6 








lOlP-A 




18.7 


1.0 


23.9 


1.0 








118P .. 




7.6 
26.4 


2.2 

1.5 


28.4 


1.1 








119P .. 






129P .. 








28.6 


1.4 








152P .. 








9.1 


1.5 








159P .. 




25.3 


2.9 


14.1 


1.6 








173P .. 




178.0 


6.8 


185.5 


3.0 








213P .. 




196.3 


5.1 


205.6 


4.7 








246P .. 




397.0 


7.4 


424.6 


7.9 








C/2005 


W2 . 










6.1 0.6 


P/2004 


Al .. 










139.1 2.0 


P/2004 


H2 .. 










6.9 1.0 


P/2004 


V5-A 


41.3 


3.1 


32.6 


1.9 






P/2004 


VRs. 






47.5 


4.6 






P/2005 


JDins 










72.6 0.7 


P/2005 


T5 .. 










2.3 0.6 


P/2005 


W3 . 










15.9 0.7 


P/2005 


Y2 .. 










162.1 1.7 



Table 6: Linear correlation coefficients and Spearman p tests used to test if color-temperature excess {T /Tbb) 
for the 11 IRS peak- up center apertures (Table Isl) is correlated with heliocentric distance, flux, and ef p (flux 
and ef p are measured at 22 or 24 pm). For the Spearman test: Z is the number of standard deviations by 
which the Spearman p deviates from its null hypothesis expected value; Pz is the probability that a value 
of Z or greater may occur with uncorrclatcd data sets; (Z) is the mean Z from Monte Carlo runs based on 
our measurements and errors. 



Parameter 


r 


Z 


Pz 


{Z) 


Th 


0.77 


2.1 


0.03 


1.8 


1 


0.09 


0.6 


0.60 


0.5 


F22 


-0.44 


1.6 


0.12 


1.5 


e/p22 


-0.44 


1.1 


0.31 


1.1 
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Appendix A. The Quantity ef p 



For observations of light scattered by comet comae, A'Hearn et al. (1984) defined the 



quantity Af p, where A is the grain albedo, and / is the areal filling factor within an observed 
aperture of radius p. The Af p parameter has been widely reported for visible light observa- 
tions of comets, and, when given in units of cm, is crudely proportional to dust production 
rate in kg s~^ ( A'Hearn et al.|[l995 ). We have defined a new parameter, e/p, that we intend 
to be an analogous quantity for observations of thermal emission from comet comae. The e in 
ef p is the effective emissivity of the dust grains; / and p remain the same as for Af p. Both 
Af p and ef p are independent of aperture size for comae with a p~^ mean surface brightness 
profile. To compute ef p from observations, we provide the formula 



efp 



hp 



(A.i: 



where I\ is the observed mean surface brightness in an aperture of radius p centered on the 
comet, Bx is the Planck function evaluated at Tc with same units as I\ (e.g., MJy sr~^), and 
Tc is the temperature of the effective continuum. In the mid-IR, Tc = 1.1 Tbb can usually 
be assumed in the absence of any color information, as discussed in §4.2[ However, it is 
important to note that efp is sensitive to the adopted T^, and variations by factors of 2 in 
efp can easily be accounted for by small variations in Tc. 

To give our SEPPCoN observations context, we present efp values based on mid-IR 
observations of several comets in the literature (Table |5|. Continuum temperatures were 
taken from the original investigations when possible; we assumed Tc = 1.1 Tbb for the 
observation of Comet 2P/Encke at 2.5 AU, and 1.08 Tbb for 103P/Hartley 2, reported by 



Sitko et al. (2010). The flux of Comet Encke's nucleus, based on a 2.3 km effective radius 



(Kelley et al. 2006), has been subtracted from the Encke photometry of Reach et al. (2007) 



and Gehrz et al. (1989) (0.04 and 1.5 Jy, respectively). For Comet Hartley 2, we computed 



a nucleus flux of 0.97 Jy assuming a 0.6 km effective radius ( Groussin et al.|[2b04 Lisse et al. 
2009). Nuclei were computed using the near-Earth asteroid thermal model of Harris (1998), 
with an IR-beaming parameter 7] = 1.0. The nucleus fluxes for the remaining comets are 
assumed to be negligible. 
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Table A. 7: efp values computed from selected comet observations in the literature, presented for comparison to the comets in our survey. The table 
has been sorted by efp. 



Comet rh X p efp a 

(AU) (iim) (km) (cm) (cm) 



Notes 



2P/Encke 2.52 23.7 18 300 
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o 



73P-B/Schwassman-Wachmann 3 

103P/Hartley 2 

73P-C/Schwassman-Wachmann 3 

2P/Encke 

C/1995 01 (Hale-Bopp) 

C/1996 B2 (Hyakutake) 

IP/Halley 

C/1995 01 (Hale-Bopp) 



Reach et al. (2007 



1.11 


11.7 


430 


744 


8 


1.06 


11.6 


330 


814 


26 


1.09 


11.7 


350 


2 518 


25 


0.38 


10.0 


2 000 


2 720 


60 


2.8 


7.8 


3 300 


52 400 


5 200 


0.34 


10.3 


14400 


53 500 


1500 


0.59 


10.3 


12 600 


127000 




0.92 


10.3 


13 400 


1280 000 


120 000 



1 week post-out 



Meech et al. 



Barker et al. 



Gehrz et al. ( M? 



(2011 



jurst, Harker et al. (2011 



(2011 



pre-perihelion, Wooden et al. ( 1999 



Mason et a l. (1998i 
Gehrz and Ney|(,1992] 



Mason et al. "(2001 



